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SUMMARY 


This paper presented analytical results from a three-dimensional finite 
element analysis, which uses a vanishing-fiber-diameter material model. One 
very attractive feature of this program is that it requires only the properties 
of the individual constituents to conduct an analysis of any metal matrix 
laminate. Additionally, PAFAC calculates the fiber stresses in each element. 
Examples are shown for specific metal matrix composites such as boron/aluminum 
and silicon-carbide/aluminum. However, the analysis and techniques described 
can be applied to any combination of matrix and continuous fiber. 

Specimen stress-strain behavior and stress at first fiber failure were 

predicted for boron/aluminum laminates containing circular holes and crack- like 

slits. The predictions compared very well with test data for laminates 

containing 0° fibers and reasonably well for [±45]2 g laminates. Mesh 

configuration was shown to have an effect on the calculation of stresses local 

to the notch. The presence of thin interface layers of matrix material had a 

significant influence on the slit tip stress state, causing sharper stress 

gradients near the notch. Interface layers reduced the slit-tip fibers stresses 

in a [±45] silicon- carbide/aluminum laminate but increased them in a [0/90] 
s s 

laminate . 

The results presented here show that the analysis of metal matrix 
composites is not trivial. Many aspects of the analysis (e.g., mesh geometry, 
interface layers, material properties) can have a significant influence on the 
stress state local to the notch. 



INTRODUCTION 


Metal matrix composites have several Inherent properties, such as high 
stiffness- to-weight ratios and high strength- to -weight ratios, which make them 
attractive for structural applications. These composites also have greater 
transverse strength, a higher operating temperature range, and better 
environmental resistance than current polymer matrix composites. Like polymer 
matrix composites, however, metal matrix composites are very notch sensitive. 
The degree of this sensitivity depends on notch size and shape, as well as the 
laminate orientation. Unlike typical polymer matrix composites, metal matrix 
composites may exhibit wide spread yielding of the matrix before laminate 
failure. To design damage- tolerant structures (or to simply understand the 
effects of fastener holes), the laminate fracture strengths must be known for a 
wide range of ply orientations, notch geometries, and loading conditions. A 
method for predicting fracture strength is needed to avoid testing all the 
laminate, notch, and loading combinations of interest. 

In metal matrix composites the matrix may yield, whereas the fibers remain 
elastic until they fracture. An analysis that models this two phase behavior 
for any specimen configuration and incorporates fiber failure was presented in 
[1]. The analysis is a three-dimensional finite element program [2] called 
PAFAC (Plastic and Failure Analysis of Composites). The PAFAC program uses the 
vanishing-fiber-diameter continuum model developed by Bahei-El-Din and Dvorak 
[3-5] to represent the essential aspects of the elastic-plastic behavior of the 
composite lamina. 

This paper will first present a brief description of the PAFAC finite 
element program and review some past work using the program. Then the results 
of some current work will be presented. The effects of mesh configuration on 
slit tip stress concentrations and fiber failure will be discussed. The effect 
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of modeling thin layers of matrix material between each orthotropic lamina will 
be discussed. 


ANALYTICAL MODELING OF FIBROUS COMPOSITES 

In metal matrix composites, the matrix yields, whereas the fibers generally 
behave elastically until they break. To model the behavior of metal matrix 
composites, the elastic -plastic behavior of the matrix and the elastic behavior 
of the fiber should be accounted for. An analysis that models this two phase 
behavior was used in the present study. This analysis was conducted with a 
three-dimensional finite element program [2] called PAFAC , which was developed 
from a program written by Bahei-El-Din et al. [ 3 , 6 ]. PAFAC uses a constant 
strain, eight-noded, hexahedral element. Each hexahiedral element represents a 
unidirectional composite material whose fibers can be oriented in the 
appropriate direction in the structural (Cartesian) coordinate system. 

Material Model 

The PAFAC program uses a vanishing- fiber-diameter material model developed 
by Bahei-El-Din and Dvorak [3-5] to represent the essential aspects of the 
elastic-plastic behavior of the composite lamina. The material model is briefly 
described in the appendix. This material model requires only the material 
properties of the fiber and matrix to predict the laminate response. The 
program is able to calculate the fiber stresses in each element from the 
material model. Further details on the material model can be found in [3-5]. 

Finite Element Model 

All the specimens analyzed contained center notches, either circular holes 
or crack-like slits. A typical specimen with a crack-like center slit is shown 
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In Figure 1. Two typical finite element meshes, a radial mesh and a rectangular 
mesh, are shown in Figures 2 and 3. Typically, a mesh consisted of 2250 nodes 
and 1392 elements, although this varied depending on the number of plies and 
specimen configuration. In all cases, only one-eighth of the specimen was 
modeled because of symmetry. A uniform stress was applied to the end of the 
specimen to simulate tensile loading. 

At least one layer of elements is needed for every ply angle in the 
laminate. For example, a [0]^ laminate could be modeled by one layer of 
elements, while a [0/±45] g would require at least three layers of elements. 
Unless otherwise stated, each ply of the laminates was modeled with one layer of 
elements . 


The following equation, developed by Johnson et al. [1], is used as the 
fiber failure criterion in the PAFAC program: 



1 


( 1 ) 


In this equation, the subscript 3 represents the fiber axial direction and 1 

represents the transverse direction. a ^ represents the current tensile 

f 

stress in the fiber and o ^ represents the current shear stress in the fiber. 
The ultimate strengths of the fiber in the tensile and shear modes are and 

<7 13 t> res P ectivel y- Previous analysis [1] did not use a mesh fine enough to 
model an element with a width of one fiber spacing. Thus, instead of using the 
ultimate tensile strength of a single fiber, an averaging technique was used in 
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[1] to determine the fiber ultimate strengths in the tensile and shear modes. 
Since two strengths are unknown, two test cases were required. The monolayer 


with a slit and a [±45]£ s specimen containing a hole were used. These two 

layups were chosen because the dominant fiber stress at the slit tip in the 

f 

monolayer is the tensile stress 0^3 > while in the [± 45 ] 2 S laminate the dominant 

f 

fiber stress is the shear stress o^. From radiographs [1], the laminate 

stress at first fiber failure could be closely estimated. The finite element 

model for each laminate configuration was loaded to the stress level where first 

f f 

fiber failure occurred. The corresponding fiber stresses a ^ and a ^ were 

recorded for the element where first fiber failure was observed. The pairs of 
f f 

o ^ and <^3 were substituted into equation ( 1 ) and the equation was solved 
for the two unknowns, ^ 3 ^ and . This gave values of 2500 MPa and 100 

MPa for ^ 3 ^ and ^ 3 ^* respectively. First fiber failure was predicted for 
boron/aluminum specimens containing either circular holes or crack-like slits in 
[ 1 ] using these values. 


NOTCH STRENGTH 

The PAFAC program was used to predict the stress -strain behavior and first 
fiber failure of boron/aluminum laminates with either circular holes or crack- 
like slits [1]. Material properties of the boron and aluminum are given in 
Table 1. This section will summarize the results presented in [1]. 

Specimens With Holes 

The overall stress-strain behavior and first fiber failure were predicted 
for [0]^, [0/145]^, [0^/145] s> boron/aluminum laminates containing 

circular holes in [1]. In Figure 4 the applied stress is plotted against the 
specimen overall strain for a [0/±45] g laminate containing a circular hole. The 
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solLd line represents the predicted response up to first fiber failure. The 
symbols represent experimental data. There is excellent correlation between the 
predicted and the experimental stress - strain behavior. First fiber failure was 
predicted to occur near laminate failure. First fiber failure was not observed 
for this laminate [1] and therefore was assumed to occur at laminate failure. 

The first fiber failure was predicted next to the hole in the 0° ply element on 
the transverse axis, as shown in the sketch. Figure 5 summarizes the first 
fiber failure predictions for the specimens containing circular holes. 
Predictions for the other layups agreed well with the experimental data. 

Specimens With Slits 

The PAFAC analysis was also used to predict first fiber failure and overall 
stress-strain behavior for several boron/aluminum laminates containing crack- 
like slits [1]. Figure 6 compares the analytical and experimental results for a 
[ 0/145 ] specimen containing a center slit. The predicted stress-strain 
response agreed very well with the experimental data. The predicted first fiber 
failure was very close to the laminate failure stress. Figure 7 summarizes the 
first fiber failure predictions for the specimens containing crack-like slits. 
Except for the [0]^, predictions for the other layups agreed well with the 
experimental data. 


EFFECT OF MESH CONFIGURATION 

It is well known that unidirectional laminates develop long, narrow plastic 
zones at the slit tip due to the intense shear stresses. These zones, parallel 
t;o the fibers, were noted by Goree and Jones [7] and are shown in Figure 8. 
However, it is not so well known that these zones of intense shear stresses 
would be present in slit laminates containing both 0° and angle plies. As shown 
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in Figure 9, Post et al. [8] found rather significant shear zones at a slit tip 
in a [0/±45] s boron/aluminum laminate. Previous PAFAC analyses [1,2] used a 
radial mesh, such as shown in Figure 2. Because of the configuration of the 
mesh above the slit tip, the radial mesh could not predict the intense shear 
stresses parallel to the fibers in the 0° direction. A rectangular mesh, with a 
rectangular slit tip, was then used [9] to capture these zones of high shear. 

In [9], five layups of silicon-carbide/aluminum ([0]g, [0^/145]^, [ 0/90 ] 2 s > 

[0/145/90] and [ ^ 1 2 s ^ tested statically to failure. The 

strain distribution ahead of the slit tip was found experimentally with a series 

of strain gages bonded ahead of the slit tip. For all layups, except the [0] o , 

o 

the yielding of the metal matrix caused the fiber stress concentration factor to 
increase with increasing load. This is contrary to behavior seen in homogeneous 
materials where yielding causes the stress concentration to drop. For the [0] o 
laminate, yielding of the matrix was found to cause a decrease in the fiber 
stress concentration. The strain distribution ahead of the slit was also 
predicted with the PAFAC analysis using a rectangular mesh such as shown in 
Figure 3, For all layups, the finite element analysis, using the rectangular 
mesh, predicted the experimental trends correctly. The predictions made with 
the radial mesh, however, were in error, predicting that the fiber stress 
concentration in the [ 0 ] g laminate would increase due to yielding. 

The effect of the mesh in the slit tip region is further investigated here. 
The two meshes shown in Figures 2 and 3, as well as the mesh in Figure 10, were 
used. The mesh shown in Figure 10 was formed to combine a curved slit tip with 
a rectangular mesh. In these figures, due to symmetry, only one quarter of the 
plan view of the mesh is shown. Each figure also shows the slit- tip detail for 
each mesh. In all three meshes, the slit- tip element has the same width 
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(dimension in the x-direction) . The slit-tip element was sized to represent one 
fiber spacing (0.192 mm). 

Notch Tip Stresses 

The three meshes were used to analyze a unidirectional boron/aluminum 
monolayer with a center slit. The material properties are given in Table 1. 

The specimen was 76 mm wide with a slit length of 19 mm. The different mesh 
configurations did not have any effect on the overall stress - strain behavior of 
specimen. The major difference was seen near the slit tip. Figure 11 compares 
the distributions of the stress in the loading direction ahead of the slit tip 
for the meshes shown in Figures 2, 3, and 10. At the slit tip, the stresses 
computed with the two rectangular meshes, for both the laminate and the fiber, 
are significantly higher than computed with the radial mesh. Stresses are 
calculated at the element centroid. Although all the slit- tip elements are the 
same size in the x-direction (parallel to the slit), the slit tip element in the 
radial mesh does not have the same dimensions in the y-direction (perpendicular 
to the slit) as in the rectangular meshes. This difference is shown in Figures 
2, 3, and 10 in the sketch of the slit tip elements. Thus the location of the 
stress for the radial mesh is slightly different than for the two rectangular 
meshes. Because of the very steep stress gradients at the slit tip, even a 
small difference in location may have a significant effect on the stress 
calculations. Away from the slit tip, all three meshes agree even though again 
the locations of the stresses are different for the radial mesh. The results 
computed using the rectangular mesh with the rectangular slit tip and with the 
curved slit tip are nearly identical. Thus, for the remainder of this paper, 
the only rectangular mesh used will be the one with the rectangular slit tip. 
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Figures 12 and 13 show the a stress contours for the radial and 

rectangular meshes shown in Figures 2 and 3, respectively. Due to symmetry 

conditions, only one quarter of the plan view of the specimen is shown in the 

stress contours. Figures 14 and 15 show the r stress contours for the 

xy 

meshes shown in Figures 2 and 3, respectively. By comparing Figure 12 to 13 and 
14 to 15, it is obvious that the radial mesh cannot model the shear band. The 
stress contours were significantly affected by the mesh configuration. The 
radial mesh (Figure 2) could not model the high shear zones present in the 
unidirectional specimen. 


Fiber Failure Strengths 

Since the rectangular mesh has a slit-tip element sized to contain one 
fiber, new values of < 7^3 and should be calculated. The averaging 

technique, as described above and in [1], was used. It was anticipated that now 
22 *" wou ld represent the ultimate tensile strength of a single fiber. Values of 
5300 MPa and 40 MPa were found for ^ 33 ^ and a 13 t ’ respectively. This value 
of: a 33 fc is rather large compared to the tensile ultimate strength of the boron 
fiber of 3160 MPa [10]. However, this value of ultimate strength for the boron 
fiber was obtained from tensile tests of unnotched 0 ° boron/aluminum specimens 
[ 10 ] in which a relatively large volume of material was subjected to a uniform 
stress. For the specimens with a notch the material next to the notch is 
subjected to a stress gradient and the peak stress exists only over a small 
volume. Thus, may exceed the fiber ultimate tensile strength. 

Using the fiber ultimate strengths of 5300 and 400 MPa in Eqn. (1) with the 
rectangular slit-tip mesh (Figure 3), first fiber failure was predicted for two 
laminates with slit-like slits previously analysed in [ 1 ] using a radial mesh. 
However, the new predictions of first fiber failure were grossly in error, 
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predicting first fiber failure at stress levels that were much too low. These 
erroneous predictions may have been due to mesh effects. The rectangular mesh 
models the shear stresses well but may not model other aspects as adequately. 
More work will be done in this area. 

EFFECT OF MATRIX INTERFACE LAYERS 

In [11] Guo, Post, and Czarnek compared Moire fringe patterns for a 

[0/±45] and a [±45/0 9 ] laminate, both center notched tensile specimens made of 
s z s 

boron/aluminum. Figure 16 compares the fringe patterns of y- displacements for 
the two specimens. The two fringe patterns have similar global deformations 
(Figure 16(a) and 16(c)). As stated in [11], this indicates that the successive 
plies of different orientations constrain each other by the interlaminar shear 
stresses acting in the matrix layers between the plies causing them to deform 
together. However, the local displacements at the slit tip (Figures 16(b) and 
16(d)) are very different. There was no evidence of delamination; therefore, 
strong shear strains must exist in the matrix layers between plies to maintain 
continuity . 

To account for this effect in the analysis, a thin isotropic layer of 
matrix material was modeled between each layer of orthotropic material. Two 
silicon-carbide/aluminum laminates were analyzed: [±45] g and [0/90]^. The 

material properties are given in Table 1. For each laminate, three cases were 
compared, no isotropic layers, one isotropic layer, and three isotropic layers 
modeling the interface layer. The model is shown schematically in Figure 17 for 
a [0/90] laminate. (The laminate in Figure 17 has only two plies for 
simplicity.) The thickness of the interface layer was chosen to be 0.03175 mm 
(typical fiber spacing for silicon-carbide/aluminum) and the fiber volume 
fraction of the orthotropic layers was adjusted to maintain a laminate fiber 
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volume fraction of 0.44. The interface shown in Figure 17 was modeled with one 
or three layers of elements. The specimen was 51-mm wide, 164-mm long with a 
slit 19-mm long in the center. 

As expected, the overall stress- strain behavior of the specimens for either 
layup were not changed by including interface layers in the model. However, the 
stress gradients in the area of the slit tip and fiber stresses in the slit tip 
element were significantly affected. 

[±45] s Laminate 

Figure 18 compares the stress contours for a • for the three models for 
the +45° layer in the [±45] g laminate at an applied stress of 60 MPa, At this 
level of applied stress there is global yielding in the specimen. Figure 18(a) 
presents the a stress contours for the [±45] laminate with no interface 

y s 

layers, 18(b) presents the contours for the laminate with one interface layers, 
and 18(c) presents contours for the laminate with three interface layers. As 
shown in the figure, adding the interface layers allowed the two layers of the 
laminate to displace relative to one another, resulting in sharper stress 
gradients above the slit tip. Figure 18 shows that not only is there a much 
sharper gradient above the slit tip but the stress contours above the slit are 
also more compressed by adding the interface layers. Similar behavior was seen 
for the shear stress components, as shown in Figure 19 for the r stress in 
the 45° layer in the [±45] g laminate. 

Including interface layers also had an effect on the fiber stress in the 
element at the slit tip. Figure 20 shows the axial component of slit- tip fiber 
stresses for the -45° layer of the [±45] g laminate at an applied stress of 60 
MPa. Here the presence of the interface layers caused the fiber stresses to 
drop significantly, about 30% for the three interface layer case. Including 
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three interface layers reduced the shear component of the fiber stress at the 
slit tip by about 3%. 

[ 0/90 ] s Laminate 

Figure 21 presents the o stress contours in the 0° layer of a [0/90] 

y s 

laminate for an applied stress of 110 MPa. A similar trend is seen here. By 
adding interface layers, the stress contours are more compressed. Also notice 
that in the area above the slit that the stresses were lowered by adding the 
interface layers. The presence of the interface layers allows the two plies to 
act more independently; thus, the 0° layer of the [0/90]^ laminate now acts more 
like a monolayer. Figure 22 shows the fiber axial stresses in the slit- tip in 
the 0° layer for the [0/90] s laminate at an applied stress of 110 MPa. For the 
[ 0/90 ] s laminate, the interface layers caused an increase of about 15% in the 
fiber stresses at the slit- tip. The interface layers allowed the 0° layers to 
experience less constraint from the 90° layers, thus producing a higher fiber 
stresses . 


CONCLUDING REMARKS 

This paper presented analytical results from a three-dimensional finite 
element analysis, which uses a vanishing-f iber-diameter material model. One 
very attractive feature of this program is that it requires only the properties 
of the individual constituents to conduct an analysis of any metal matrix 
laminate. Additionally, PAFAC calculates the fiber stresses In each element. 
Examples are shown for specific metal matrix composites such as boron/aluminum 
and s ilicon- carbide/aluminum . However, the analysis and techniques described 
can be applied to any combination of matrix and continuous fiber. 



Specimen stress-strain behavior and stress at first fiber failure were 
predicted for boron/aluminum laminates containing circular holes and crack- like 
slits. The predicted stress- strain response and stress at first fiber failure 
compared very well with test data from laminates containing 0° fibers and 
reasonably well for [±45]2 S laminates. 

It was shown that the mesh configuration can have an effect on the 
calculation of the stress distribution in the notch region. The kinematics of 
the problem must be considered when determining the mesh configuration. 

The presence of thin interface layers of matrix material had a significant 

influence on the calculated slit tip stress state in silicon-carbide/aluminum 

laminates F causing sharper stress gradients and allowing each ply to act more 

independently. Interface layers reduced the slit- tip fibers stresses by 30% in 

a [±45] s silicon-carbide/aluminum laminate but increased the slit- tip fiber 

stress by 15% in a [0/90] silicon-carbide/aluminum laminate. 

s 

The results presented here show that the analysis of metal matrix 
composites is by no means trivial. Many aspects of the analysis (e.g., mesh 
geometry, interface layers, material properties) can have a significant 
influence on the stress state local to the notch. 
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APPENDIX - MATERIAL MODEL OF PAFAC PROGRAM 


The PAFAC program uses a vanishing- f iber -diameter material model developed 
by Bahei-El-Din and Dvorak [3-5] to represent the essential aspects of the 
elastic-plastic behavior of the composite lamina. The material model is briefly 
described in this appendix. 

The model consists of an elastic-plastic matrix unidirectionally reinforced 
by continuous elastic fibers. Both constituents are assumed to be homogeneous 
and isotropic. The fibers are assumed to have a very small diameter, so that 
although the fibers occupy a finite volume fraction of the composite, they do 
not interfere with matrix deformation in the two transverse directions, but only 
in the axial (fiber) direction. This assumption is generally referred to as the 
vanishing-fiber-diameter model. Figure 22 shows a schematic of this lamina 
model. In Figure 22, the x^, x^ , and x^ axes are the element coordinate system. 
In Figure 22, the fiber (axial) direction is parallel to the x^-axis, and the 
x^- and x^-axes represent the transverse directions. The lamina model can also 
be; represented by parallel fiber and matrix bars or plates with axial coupling, 
as illustrated in Figure 22. 

If the Cartesian coordinates are chosen so the x^ coincides with the fiber 
direction, the second-order tensors of the independent stress and strain 
components, a and e, are expressed as 

° “ [a ll °22 a 33 °\2 a 13 CT 23 j 

6 = le ll € 22 € 33 7 12 7 13 7 23 1 

where 7 ^ -2c_ (i, j * 1, 2, 3; i*j) are the engineering shear strain 
components . 
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For equilibrium and compatibility, several requirements are imposed on the 
material model shown in Figure 22. The only constraint in the model is in the 
axial (fiber) direction; the matrix and fiber must deform equally. Thus, 

e 33 ” **33^ _ ^33^ 

The stress average in each constituent can be related to the overall 
composite stress o in the axial (fiber) direction as follows: 

"33 " v f (a 33 ) f + v m <a 33*in 

A bar over a symbol indicates overall composite stress or strain, and the 

subscripts f and m denote quantities related to the fiber and matrix. The 

volume fractions v_ and v are such that v~ + v =1. In this material 

r m f m 

model, the other stress components in each constituent are assumed to be uniform 
and to obey the following equilibrium equations: 

"ll = ^l^f = ^l^m 
"22 " (a 22*f “ (a 22 ) m 
"12 - ( "l2>f - ^"l2 ) m 
"13 = ( "i3>f = < a i3>« 

"23 “ ( "23 ) f “ ^ a 23^m 

The other strain components can be related to the overall strain, i as 
follows: 


ij 


• • ) c + v ( £ • • ) 

f 1J f m ij m 


(ij * 33) 
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Since the fibers are elastic up to failure, the inelastic strains of the 
lamina are caused by matrix deformation. Because the fiber imposes an elastic 
constraint on the matrix which affects the shape of the lamina yield surface, 
additional kinematic components appear in the hardening rule of the lamina and 
influence the magnitude of the overall plastic strains. All aspects of the 
yield behavior were examined and accounted for in the formulation of the lamina 
constitutive equations. These equations are explicitly described in [2,3]. The 
stress -strain curve of the matrix material was modeled with a Ramberg-Osgood 
equation. The slope of the effective stress -effective plastic strain curve H 
was calculated using the following equation: 


H - (3 G/n/ic) (® ys /% ff ) <n ' 1> 


m 


where a e ff is the effective stress in the matrix material, <7 is the yield 
stress of the matrix material, G is the shear modulus of the matrix material, 
and n and k are determined by fitting experimental data. For the 6061 


aluminum matrix material, G - 27200 MPa, a 
0.05. 


ys 


34 . 5 MPa , n - 5 . 5 , and k 
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Table 1 - 

Constituent Properties 


Material 

Modulus , GPa 

Poisson's Ratio 

Yield Strength, MPa 

Aluminum, 6061 

72.3 

0.33 

34.5 

Boron 

400 

0.13 


Silicon- Carbide 

340 

0.25 







Figure 1. Typical specimen configuration with a crack-like slit. 




Figure 2. Radial finite element mesh. 2a - 19 mm, W - 78 mm, and H - 101 
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Figure 4. Stress-strain response to failure of [0/±45] g laminate with circular 
hole. 2a - 25.4 mm; w - 102 mm. 
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® Stress strain response to failure of [0/±45] s laminate with crack- 

like slit. 2a - 10 mm; w - 102 mm. 
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compared with laminate ultimate stress for 


High shear stress bands 
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Figure 8. Brittle lacquer crack pattern in boron/aluminum monolayer at 
12,5000 N. Goree and Jones [7]. 
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Figure 10. Rectangular mesh with curved slit tip. 2a - 19 mm, W - 78 mm, and H 
- 101 mm. 
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Figure 11. Fiber and laminate stress distributions ahead of slit tip for 
different meshes. 


<D 

_cd 

o 

o 







II 

cd 

02 


II 

£ 


cd 

n, 

2 

c 

d> 

cd 

co 

<u 

CO 

7} 

CD 

-t-> 

71 

r-H 

< 




^ CO 00 o 

O O O 'TH 



30 


Figure 12. a stress contours for boron/aluminum monolayer using radial mesh. 
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Figure 13. a y stress contours for boron/aluminum monolayer using rectangular 
mesh. 



Figure 14. r stress contours for boron/aluminum monolayer using radial mesh. 
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Figure 15. r xy s ^ ess contours for boron/ aluminum monolayer using rectangular 



(a) Global pattern, [0/±45] 


(b) Local pattern at slit tip, [0/±45] 



(e) Global pattern, [±45/0 2 ] s . (d) Local pattern at slit tip] [±45/0 ? ] 


I’if’urc' 1ft. Moire fringe patterns for a [0/±45] g and a [±45/0 2 ] g boron/aluminum 
tensile specimen with a center slit. 2a = 5.5 mm, W = 14.7 mm. 

Guo, Post, and Czarnek [11]. 
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Figure 17. Interface 
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Figure 23. Element material model. 
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